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Degenerate optomechanical parametric oscillators are optical resonators in which a mechanical 
degree of freedom is coupled to a cavity mode that is nonlinearly amplified via parametric down- 
conversion of an external pumping laser. Below a critical pumping power the down-converted field 
is purely quantum-mechanical, making the theoretical description of such systems very challenging. 

Here we introduce a theoretical approach that is capable of describing this regime, even at the 
critical point itself. We find that the down-converted field can induce significant mechanical cooling 
and identify the process responsible of this as a cooling-by-heating mechanism. Moreover, we show 
that, contrary to naive expectations and semi-classical predictions, cooling is not optimal at the 
critical point, where the photon number is largest. Our approach opens the possibility for analyzing 
further hybrid dissipative quantum systems in the vicinity of critical points. 

PACS numbers: 42.65.Yj,42.50.Wk,42.50.Lc,03.65.Yz 


I. INTRODUCTION 

Degenerate optical parametric oscillators (DOPOs) 
consist of a driven optical cavity containing a crystal with 
second-order optical nonlinearity nHa. Down-conversion 
in the crystal can generate a field at half the frequency 
of the driving laser and classical electrodynamics pre¬ 
dicts that such field will start oscillating inside the cav¬ 
ity only if the external laser power exceeds some thresh¬ 
old value, where the nonlinear gain can compensate for 
the cavity losses. A fully quantum-mechanical theory, 
on the other hand, reveals that even below threshold the 
down-converted field is not vacuum, but a squeezed field 
whose quantum correlations increase as the threshold is 
approached. 

Recent developments in the fabrication of crystalline 
whispering gallery mode (cWGM) resonators [5HT7] have 
opened the way to study the intracavity interplay be¬ 
tween down-conversion and optomechanics m, a setup 
that we will refer to as degenerate optomeehanieal para- 
metrie oseillator (DOMPO). So far, it has been shown 
that the presence of down-conversion in an optomechan¬ 
ical cavity can help enhancing mechanical cooling m, 
normal mode splitting [20], sensitivity in position mea¬ 
surements [21], or even bringing optomechanics close to 
the strong coupling regime with additional bath engi¬ 
neering [22]. In all these works, however, the nonlinear 
crystal is operated as a parametric amplifier, providing 
a nonlinear gain to some external field that is injected 
in the cavity at the down-converted frequency {stimu¬ 
lated down-conversion). In contrast, the description of 
the interaction between the field generated via sponta¬ 
neous down-conversion and the mechanical mode is much 
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more challenging, since (below threshold) the former is 
purely quantum mechanical [23] , so that the optomechan¬ 
ical coupling cannot be linearized and does not admit a 
simple Gaussian description. 

In this work we provide a theory for the DOMPO 
which can be trusted all the way to threshold, and is 
obtained by combining traditional adiabatic elimination 
techniques with our recently developed self-consistent 
Mori projector (c-MoP) theory [24l [25]. To this end 
we first introduce the master equation which models the 
DOMPO, and perform an adiabatic elimination of the 
optical modes by neglecting the mechanical backaction. 
The mechanical state is found to stay approximately 
thermal for parameters compatible with current cWGM 
resonators, with an effective temperature dependent on 
the steady-state value and two-time correlation function 
of the down-converted photon number, which we derive 
in two ways. First by treating the pump mode as a clas¬ 
sical field (semi-classical approach), allowing us to obtain 
simple analytical expressions and provide a physical ex¬ 
planation for the regions of significant cooling, showing 
that the system provides a realistic implementation of the 
cooling-by-heating mechanism [26] below threshold. Sec¬ 
ond, by using c-MoP theory on the optical dynamics to 
find reliable results at threshold and justify the absence 
of the mechanical backaction onto the optics. Remark¬ 
ably, this accurate approach allows us to prove that the 
semi-classical predictions break down when working very 
close to threshold, where cooling is shown to disappear. 
These results might have strong implications not only for 
future analysis, but also for previous results which make 
use of semi-classical approaches while working very close 
to threshold [22]. In the final section we apply c-MoP 
theory on the full optomechanical problem, and identify 
the region of the parameter space where the mechanical 
backaction on the optics is negligible (which contains the 
parameters of our interest). 
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II. THE DEGENERATE OPTOMECHANICAL 
PARAMETRIC OSCILLATOR. 

The system we consider can be schematically repre¬ 
sented as in Fig. A crystal with second-order optical 
nonlinearity is shared by two cavities with relevant reso¬ 
nances at frequencies uOp (pump) and ujg ^ uOpj2 {signal). 
The pump cavity is driven by a resonant laser, so that 
photons in the signal cavity can be generated via spon¬ 
taneous down-conversion mill. In addition, one of the 
mirrors of the signal cavity can oscillate, and is therefore 
optomechanically coupled to the down-converted field via 
radiation pressure [iH]. Let us define annihilation oper¬ 
ators {cij}j=p,s,m for the pump (p), signal (s), and me¬ 
chanical (m) modes. Including losses of the optical modes 
at rate 70 (assumed the same for pump and signal with¬ 
out loss of generality), as well as the irreversible energy 
exchange of the mechanical mode with its thermal envi¬ 
ronment at rate 7 ^ (with which, in the absence of light, it 
is in thermal equilibrium with nth phonons), the master 
equation governing the evolution of the DOMPO’s state 
p can be written as 

P = -'AH, p] + loHa, [p] + 'lo'Da^ [p] ( 1 ) 

"t“ Tm(^th "t“ VHam [/?]• 

We have defined Lindblad superoperators Vj[-] = 

2J(-)J”*" — J(') — (•)J, and the Hamiltonian 

H — -Hopt “ 1 “ ^m^\n^Ta (^m ^m) •> ( 2 ) 

where we normalize the optomechanical coupling 770 M 
to the frequency of the mechanical oscillation The 
optical Hamiltonian can be written in a picture rotating 
at the laser frequency as 

L^opt — OpOg) , (3) 

where = ujs — ujp12 is the detuning of the signal 
mode (which we will take positive in this work), x/2 is 
the down-conversion rate, and Cp is proportional to the 
square root of the injected laser power. 

In the classical limit, the steady-state phase diagram 
of the DOMPO features a variety of phases [23]. Here 
we focus on the regime where the state of the signal 
field is fully quantum, i.e. where the trivial solution 
{os) = 0 is the only stable one, henceforth referred to 
as the monostable phase^ which requires two conditions. 
First, defining the injection parameter a = Cpxl^fi and 
the normalized detuning A = A 5 / 70 , the trivial solution 
becomes unstable in favour of a nontrivial one {os) 7 ^ 0 
for a > Vl -h A^ [ 23 ]. Hence, we write a = Vl + A^ x 
and focus on the x G [0, 1] region. The second condi¬ 
tion, 4ClAr]Qy^/r]^Q < 1, guarantees that the nontrivial 
solution does not enter the x G [0, 1] region [23]. In this 
expression we have introduced the dimensionless down- 
conversion coupling ?7 dc = x/lo as well as the sideband- 
resolution parameter Q = Clrn/xo- 
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FIG. 1. (Color online) Sketch of the degenerate optome¬ 
chanical parametric oscillator. We take crystalline whisper¬ 
ing gallery mode resonators ISH 13 as a reference for our 
choice of parameters, but electromechanical implementations 
are also possible using superconducting circuits in the de¬ 
generate parametric oscillation configuration m coupled to 
drum-shaped oscillating capacitors |28ti3Q] (see also El). 


We emphasize that the vanishing signal field amplitude 
excludes the possibility of using a linearization approach 
similar to those applied in [TilllllH31|. In the follow¬ 
ing, we provide a theory that works in all the x G [0,1] 
region, and use it to predict the action of the down- 
converted field on the mechanical state. 


HI. EFFECTIVE MECHANICAL DYNAMICS 


Despite the complexity of the problem, we remarkably 
find with the help of c-MoP theory that for typical system 
parameters the optical modes do not receive considerable 
backaction from the mechanics. This property, which 


we justify in Sec. [VT| allows us to simplify the problem 
significantly via an adiabatic elimination of the optical 
modes [32] I37II4Q] leading to an effective master equation 
for the reduced mechanical state Pm(^)- As we show in 
Appendix the mechanical steady state can then be 
approximated by a thermal state (displaced by (a^) « 
VoM^s) characterized by its phonon number 


^th + I + r>i ^ 

i + (r_-r+) ^ F 


^FL , 


( 4 ) 


where F = r_— F^ is the cooling efficiency and 
npL = r+/r the fundamental limit for the phonon 
number. All the information about the optical modes 
is contained in the heating and cooling rates T± = 
CRejyo dr exp(^'if7^r) 5(r)} through the optical 
correlation function 


s(r) = tr{4as e^°p‘^a7s Popt} - (5) 

where Ns = tr{a|asPopt} is the signal photon num¬ 
ber and C = cooperativity with Q = 

O^m/lm fhe mechanical quality factor. Here, Popt is the 
steady state of the optical Liouvillian 

>^opt[’] = ~i[L^optG] A ^ ^ 7oL^aj[’]: (b) 

j=s,p 
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FIG. 2. (Color online) Steady-state phonon number as a 
function of A and x, as obtained from a semi-classical descrip¬ 
tion of the DOPO and a thermal phonon number nth = 100 
(achievable either with cryogenics [41] or via a standard side¬ 
band pre-cooling stage). We choose typical parameters for 
cWGM resonators [TOj: Q = 10, ?7om = 10“^, and Q = 10®. 
Note that typical optical decay rates are on the MHz range, al¬ 
though only normalized quantites are relevant for our results. 
The signal steady-state photon numbers Ns{x) corresponding 
to the ticked x-values are shown in the upper axis, showing 
that cooling is effective even with just ~ 100 photons. On 
the other hand, c-MoP theory has allowed us to prove that 
the cooling region closer to threshold disappears once such a 
more accurate approach is used, see Fig. [^ 


that is, £opt[/0opt] = 0. 


In the following we study the behaviour of the steady- 
state phonon number as we approach the DOMPO’s 
threshold. From Eq. it is clear that optimal cool¬ 
ing is then found by simultaneously maximizing F and 
minimizing npL- 

The nonlinear nature of the parametric down- 
conversion process in Eq. and a potential back action 
of the mechanical mode preclude an exact treatment of 
the optical correlation function in Eq. (§. To get sim¬ 
ple analytic expressions that enable physical insight, we 
first apply standard linearization to the optical problem, 
which we denote by semi-classical approach and has been 
the method of choice in previous works [T9]-[22]. Next, 
applying c-MoP theory [25] we show the failure of the 
semi-classical approach close to the critical point and find 
more accurate expressions at criticality, which will also 
allow us to justify the adiabatic elimination of the optical 
modes. 


IV. SEMI-CLASSICAL APPROACH 


Below threshold, the linearization of the DOPO is ac¬ 
complished by treating the pump mode as a classical sta¬ 
tionary source, that is, by performing the replacement 
ap Cp/^Q [IHl]. Within this approximation, the optical 
problem is governed by a Gaussian single-mode Liouvil- 
lian 


7o ^^optH = -i[l^alas+a{al'^ -al)/2,-]+Va,[], (7) 


from which any correlation function can be easily found, 
allowing us to obtain analytical expressions for the rel¬ 
evant quantities in Eq. as we show in Appendix]^ 
Eor the fundamental limit, we find 

npL = [4 + (n-2A)2]/8nA, (8) 

while the cooling efficiency can be written as 

T = (9) 


where we have defined the function 




(4 + Fi2)[^4 + 16(1 + ^2^)2 + 8112(1 _ ^2^)] ’ 


and a parameter (5eff = \//S? — that will be shown later 
to play the important role of an effective optical detuning. 
The photon number Ns{x) = x^/(2 — 2x^) is fully due to 
quantum fluctuations and increases hyperbolically until 
threshold x = 1 where it diverges in this semi-classical 
approach. 

In Eig.j^we show the steady-state phonon number as a 
function of the two control parameters, detuning A and 
distance to threshold x, fixing the rest of parameters to 
typical values of cWGM resonators laiio]. There are two 
regions where significant cooling effects appear. One is 
in the vicinity of the threshold point and can be traced 
back to the vast increase of the photon number Ng which 
makes F ^ 1 for virtually any value of the remaining 
parameters. However, as we will show below with the 
c-MoP approach, so close to threshold this semi-classical 
approach breaks down, hence rendering this prediction 
incorrect. 

The other region, which turns out to be of major signif¬ 
icance when aiming for optimal cooling, corresponds to 
(5eff ~ ^4/2, see the black solid line in Fig. The c-MoP 
approach will confirm this prediction in the next section. 
Moreover, it can be understood in physical terms by mov¬ 
ing to a new picture defined by the squeezing operator 
S{r) = exp[—ir(a|^ +a^)/2] with tanh2r = a/A (note 
that such transformation requires A > cr, which corre¬ 
sponds in Fig. to the region above the black dashed 
line). This transformation diagonalizes the Hamiltonian 
in the optical Liouvillian 0 . so that, defining the pa¬ 
rameters TVeff = (A/(5eff “ l)/2 and M = cr/2^eff5 tfie 
transformed state p = S'^{r)pS{r) evolves according to 

7o ^dtp = [—ic^effujas, p] + (A^eff + [p] + [p] 

-[inal^a^,p] + {Q./Q){{nth + l)X>a^[p] + [p]} 

+ [Q7]omM{ ala'll - afam),p], (10) 
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within the rotating-wave approximation, valid under the 
conditions ^ a and t^omA <C see Appendix 
Therefore, in this picture the signal field is turned into 
a bosonic mode with oscillation frequency (5eff and ther¬ 
mal occupation TVeff- The optomechanical coupling is 
dressed by the squeezing parameter M, similarly to the 
dressing by the intracavity field amplitude in standard 
sideband cooling [32 SO]. However, at difference with 
that case the interaction exchanges phonons with pairs 
of photons (rather than single photons), thus explaining 
why (5eff = f^/2 is the resonance condition for cooling. 
Within this condition, assuming 27Veff ^ 1 and ^ 4, 
we furthermore find T ^ 2CM^7Veff- The cooling effi¬ 
ciency r thus receives an additional contribution 27Veff 
from the effective thermal photon number, which is a di¬ 
rect consequence of the nonlinear nature of the effective 
optomechanical coupling in (10) that cannot be found in 
standard sideband cooling, as we discuss in Appendix 
0 This represents a natural example of the so-called 
cooling-by-heating effect [26], where heating up the opti¬ 
cal field can contribute to making optomechanical cooling 
more efficient. However, as it is well known from stan¬ 
dard sideband cooling [42], thermal photons also con¬ 
tribute to the fundamental limit, which indeed can be 
approximated by npL ~ Aeff/2 in our scenario. When 
the term rith/T dominates over npL in Eq. (|^ the ther¬ 
mal optical background TVeff can then be interpreted as 
“good noise”, while as soon as the fundamental limit is 
reached it becomes “bad noise” and heats up the me¬ 
chanical motion (see the thick black curve in Fig. |^. We 
discuss this phenomenon in more detail in Appendix [C] 
It is interesting to examine the limits that this “bad 
noise” imposes on cooling. For the parameters of Fig. [^ 
compatible with current cWGM resonators, it is found 
that the cooling efficiency and the fundamental limit can 
be simultaneously optimized to about F « 85 and npL ~ 
3 as a function of the experimentally-tunable parameters 
A and a, precluding ground-state cooling. However, the 
foreseeable increase of the mechanical quality factor and 
the optomechanical coupling by one order of magnitude 
in next-generation cWGM resonators m will improve 
these numbers to F 400 and npL ~ 0.15. 


V. C-MOP APPROACH 

The semi-classical approach has allowed us to get an¬ 
alytical and physical insight into the problem. It is how¬ 
ever well known that this approximation fails close to 
the critical point, although there is no systematic way of 
checking where exactly within the semi-classical formal¬ 
ism itself. Hence, to determine where it exactly breaks 
down and to find more accurate results for those param¬ 
eters, we make use of the recently developed c-MoP tech¬ 
nique [24l |25] , which allows finding reduced equations for 
the constituent parts of a composite system, even in sit¬ 
uations where there is significant backaction among its 
parts and no time-scale separation between their dynam- 
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FIG. 3. (Color online) Terms nth/F and npL contribut¬ 
ing to the steady-state phonon number 0 as a function of 
the distance-to-threshold parameter x. We fix the detun¬ 
ing to A = 75, corresponding to the straight thin solid line 
in Fig. [^ taking gr>c = 0.01 [9] and the rest of parame¬ 
ters as in that figure. The inset displays a close up of the 
steady-state phonon number fim for x G [0.999,1], but with¬ 
out applying the rotating-wave and Markov approximations, 
see |A5| Note that c-MoP gives finite results equivalent to 
those found within the rotating-wave approximation, while 
the semi-classical predictions diverge at threshold, which can 
be taken as further evidence that the theory breaks down 
there. 


ics. 

For parameters compatible with cWGM resonators, 
the theory is already regularized by using c-MoP only 
in the optical problem (DOPO), which provides a more 
accurate description for the optical correlation function 
(§ and photon number that enter the effective mechan¬ 
ical dynamics. The application of c-MoP to the DOPO 
has been detailed in [25] , but we review its most relevant 
steps for completeness in Appendix[Dj Specifically, we use 
a combination of c-MoP and a Gaussian-state approx¬ 
imation, which provides an efficient and accurate tool 
capable of regularizing the divergencies and unphysical 
predictions of the semi-classical approach. In particular, 
we show in Appendix that at threshold the decay rate 
of the optical correlator scales as 7opt oc 7 o^dc( 1 + A), 
while the photon number as oc (1 + A)/? 7 dc 5 in con¬ 
trast to semi-classical results in which the former goes to 
zero while the latter diverges. 

We show a very representative case for the phonon 
number as a function of the distance to threshold 
X in Fig. [^ The method finds perfect agreement with 
the semi-classical predictions sufficiently below thresh¬ 
old, in particular verifying the cooling-by-heating effect 
presented above. Most importantly, we find that the ab¬ 
solute minimum phonon number is indeed reached when 
the resonance condition Jeff = is met. On the other 
hand, close to threshold we find a significant correction 
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to the semi-classical predictions for the fundamental limit 
npL- In particular, while this is independent of the dis¬ 
tance to threshold x in the semi-classical picture, c-MoP 
shows that it actually increases very rapidly as the crit¬ 
ical point is approached, and hence no cooling is found 
no matter how much the efficiency P is increased. This 
is consistent with the fact that when A < a (as hap¬ 
pens at threshold), Jeff becomes imaginary and there is 
no resonance for the optomechanical interaction. 


VI. ABSENCE OF MECHANICAL 
BACKACTION ONTO THE OPTICS 


The adiabatic elimination of the optical fields which we 
have used throughout the work relies on the time-scale 
separation between the optical and mechanical degrees 
of freedom. In particular, such an approach neglects me¬ 
chanical backaction onto the optics, which is a good ap¬ 
proximation as long as the rate of any mechanical pertur¬ 
bation is much smaller than the intrinsic relaxation rate 
of the optics 7 opt- Far from the critical point the optical 
relaxation rate is 70 , which usually dominates over any 
other rate in the system. However, as the critical point is 
approached the DOPO dynamics exhibits a critical slow¬ 
ing down, and its relaxation rate becomes smaller and 
smaller. Hence, in our work that considers parameters 
close to threshold, it is very important to check that the 
desired time-scale separation is present. 


An intuitive argument supporting such a time-scale 
separation follows from relating the mechanical back- 
action rate with the optical frequency-shift induced by 
the optomechanical interaction, 7 ba ck = 70 M^m(^m) = 
2r]QM^m^s^ where we have used (A 8 ). Hence, using 
the scaling of Ng and 7 opt obtained in the previous sec¬ 
tion at threshold, the condition 7 back ^ 7 opt becomes 
20 (? 7 om/ 7 dc)^ ^ 1 ? which is very well satisfied for the 
parameters that we work with. Moreover, note that this 
condition is automatically satisfied when working within 
the monostability condition 4Af](7oM/7Dc)^ < 1 as long 
as A ^ 1 / 2 . 


We can set more rigorous bounds to the region where 
mechanical backaction is negligible by using c-MoP the¬ 
ory [24l |25], since, in contrast to adiabatic elimination 
methods, it does not rely on the concept of time-scale 
separation or absence of backaction effects. Hence, we 
apply this theory to the DOMPO system by using the 
time-dependent self-consistent Mori projectors = 

Popt{t)<S)tTopt{-} and = tr^{-}( 8 )pm( 0 , that is, using 
a bipartition “optics (8) mechanics” for the system. This 
approach will allow us to identify the terms contributing 
to the mechanical backaction and find upper bounds to 
their scaling. Following the procedure introduced in pre¬ 
vious works [24l [25] , the c-MoP equations for the reduced 
optical and mechanical states in the asymptotic t ^ 00 


or steady-state limit are easily found to read 


dpn 


dt 


— 0 — ^mPm T ) [^m: Pm] ( 11 a) 


- 


Pom 


poo 

„/ dTe‘^'^'"{5XmPmS{T) -Yi.c.) 

JO 


dpo\)t 


dt 

~ ^mPOM 


— 0 — >^optPopt T l^mPOM (^m) 5 Popt] (Hb) 


pOO 

/ dre^°p‘^((5nsPoptSm('r) - H.C.) 
Jo 


where, for any operator A, we have introduced the usual 
definitions {A) = tr{A(ps(g)pp)} and SA = A — (A), while 
6ns = ~ s{r) is the usual optical correlation 

function ([^, and we have defined the mechanical corre¬ 
lation function 


^ ^^mPm} — ^ 


_ J-inm-7m)r 


( 12 ) 


+ 




{{5al^Sa„^) + {Sa^m)) + H.c. 


where the final expression is easily found by following a 
similar procedure as the one shown in Appendix for 
the semi-classical optical correlation function, since C^n 
is quadratic. 


The last two terms on the right-hand side of Eq. (11b) 
account for the mechanical backaction on the optics, 
and we proceed now to bound their effect. The sec¬ 
ond to last term describes precisely the optical detuning 
^mVouixm) induced by the optomechanical interaction 
that we already discussed at the beginning of the sec¬ 
tion, concluding that it is negligible within the parame¬ 
ter regime we work with. Then we focus on the last term 
in Eq. (11b), the Born term. Eor this purpose we de¬ 
rive the steady-state equation for the moment (a|^) (note 
that the equation of motion of (a|as) receives no explicit 
mechanical backaction even within c-MoP theory) which 
reads 


0 — fl’opt{<^s >^optPopt} (^m) ) (I^) 

poo 

+ ‘J^mVoM / t^'rRe{s„(T)}tropt {ai^e^°p*^[a7s>Popt]} 
Jo 

pOO 

-{-2iCtmVoM I I^{'^m('^)} fl'opt ^ {J^s^Popt}}? 

Jo 


where {*, •} denotes the anticommutator. Note first that 
the correlation functions tropt Popt]} and 

f^opt Popt}}: which are of similar struc¬ 

ture as the optical correlation function s(r), decay to 
zero at a rate 7 opt- Next, we derive upper bounds for the 
last two terms in Eq. (13). Eor the second to last term 
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we find 

2^1 


al e 


ala. 


pOC 

VoM / drRe{sm{r)}tTopt{ 

Jo 

on2 2 tl'opt {<^1 Popt] } 

^mVOM 


j Popt] } 


< 


Tm H“ Topt 

^ ^ ^ 

"^back 


Tm H“ Topt 


(14) 


where in the last step we have used s^(0) = 1 + 
{6al^6am) + (^^m) ^ {dal^dam) = Pm (note that we 
expect the mechanical state to stay approximately ther¬ 
mal, and hence {Sa‘^) « 0) and tropt {n|^[n|ns,Popt]} = 
([a|^, a|ag]) = 2(a|^). Similarly, for the last term in Eq. 
(13) we find 


practice backaction should be negligible even in a much 
broader region of the parameter space. 

Overall, c-MoP theory has allowed us to quantify the 
mechanical backaction on the optics in a rigorous man¬ 
ner. We have obtained very conservative bounds that the 
system parameters must satisfy in order for such backac¬ 
tion to be negligible, showing that this is indeed the case 
for the parameters used in our work, which are compati¬ 
ble with an implementation based on cWGM resonators. 
It is however foreseeable that such devices, as well as their 
electromechanical counterparts, will be able to study re¬ 
gions where backaction is significant, in which case the 
c-MoP approach presented in this section will be very 
useful. 


VII. CONCLUSIONS 


poo 

|2iU^7^^M / drlm{srn{r)}tTopt{cil‘^e^"^"'^{dns,popt}} 
Jo 


< 


^2 ^2 i'^opt {<^1 ^ {^^s^Popt}} 


2^mPOM" 


Tm “1“ Topt 


< |tropt {aJ^{(5n5,popt}}| 

7m + 7opt 

Tm Topt 


(15) 


^back 


where for the last expression we have used 
tropt_{4q<5n5,Popt}} = 2((afa^)-(at2)(atas)+(aJ2)) « 
2(2Vs_f l)(a|^), within the Gaussian state approxima¬ 
tion ( |B6| ), that is, (a|^as) ~ 3(a|^)(a|as), noting that 
{as) =0 below threshold. 

A sufficient condition for mechanical backaction to 
be negligible is then Tback^Tback < Topt- We pass to 
check whether this is the case in our work. Note first 
that 7 ^ ^ 7opt even at threshold, since 7opt/7m ^ 
7oPdc(1 + A)/7m ^ 1 for the parameters we are in¬ 
terested in. Using the scalings 7opt oc 7 oPdc( 1 + A) and 
Ns oc (1 + A)/7 dc at threshold (where these bounds 
are the tightest), we can then write the conditions under 
which backaction is negligible as 


7Lck _ . 

7opt ^Dc(l + ^)^ 
7back ^ ^ 

7opt ^Dc(l + ^) 


(16a) 

(16b) 


For the parameter set of Fig. 2 these lead to the condi¬ 
tions fim ^ 100(1 + A)^ and 1 + A ^ 1, respectively. 
For the large values of A that we use during most of 
the work, these conditions are very well satisfied. For 
small A they seem to be too tight, but we need to stress 
here that we have been extremely conservative when es¬ 
timating the Born terms (14) and (15), meaning that in 


By exploiting adiabatic elimination techniques, semi- 
classical methods, and c-MoP theory, we have provided 
a theoretical analysis of the DOMPO which works even 
at the critical point. We have focused on the region 
where the optical field is fully quantum, showing that 
such a quantum-correlated field with no coherent compo¬ 
nent can induce significant mechanical cooling through a 
cooling-by-heating mechanism. c-MoP techniques have 
allowed us to check the validity of the optical adiabatic 
elimination as well as the semi-classical approximation, 
whose predictions have indeed been shown to break down 
at threshold, showing the potential of c-MoP to treat 
dissipative quantum-optical problems in the vicinity of 
critical points. 
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Appendix A: Elimination of the optical modes 

Here we present a derivation of the effective mechan¬ 
ical master equation leading to the phonon number of 
Eq. (4) in the main text. Our starting point is the mas¬ 
ter equation governing the evolution of the state p{t) of 
the DOMPO, Eq. (1), which for convenience we rewrite 
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here as 


dtP — ^opt[p] + ^m[p] + ^om[p] (A1) 


with 


^opt [p] 

jCm[p] 


= -iA^aJa* + - H.c.j 

+ 7o2^ap [p] + 7o2^a3 W, 

= [ p\ 

“i“ 'Tm('^th “1“ [P] [p] 5 


>Com[p] = +a]^),/)]. 


^P 

(A2a) 

(A2b) 

(A2c) 


All the quantities are defined in the main text, and we 
remind the notation Vjlp] = 2JpJ^ — Jp — pJ^J for 
superoperators in Lindblad form. 

In order to eliminate the optical modes and find an 
effective master equation for the mechanical state pm{t) 
we proceed as follows. We first define the projector su¬ 
peroperator V['] = Popt^fropt{'} whose action on the full 
state p{t) of the DOMPO is V[p{t)] = Popt^pm{t). Here, 
Popt is the steady state of the optical Liouvillian, that is, 
£opt[Popt] = 0- Applying this superoperator and its com¬ 
plement 1 — V onto the master equation, and formally 
integrating the latter, we obtain an exact equation of 
motion for pm{t)^ the so-called Nakajima-Zwanzig equa¬ 
tion miisH]. Such equation is not solvable, and therefore 
we apply a Born approximation which takes into account 
terms up to second order in the optomechanical interac¬ 
tion. The resulting equation reads 


Pm(^) — “1“ ^^mPOM^s Pm{t^] 


(A3) 


^m^OM 


•f 


dr e‘^^'"[xmpm{t - t)s(t) - H.C.] 


mechanical moments form a closed linear algebraic set 
0= ( if^rn 'Tm)(^m) "T ^^mPOM-^s (A5a) 

-^mVoMRe{do}(Xm), 

0= 7m(nth - (SalnSam)) (A5b) 

-^mVOMRe{(d+ -d_){Sal^5a^) + {d*_ - d+){Sal^) - }, 

0 = {-Rim - lm){ 5 a^m) (^Sc) 

- dX){5al^5am) + (d+ - d)l)((5a^) + d_], 

where we used the abbreviations {A) = tr{Ap^}, 5A = 
A — (A), and 

POO 

dp = drs{r)^ (A6a) 

Jo 

pOO 

d±= dr s{t). (A6b) 

Jo 


These equations can be solved for the steady-state mo¬ 
ments as functions of the optical photon number Ng and 
correlation function s{r) without the need of further ap¬ 
proximations. However, in order to obtain more physical 
insight into the mechanical steady state pm we apply both 
Markov approximation and a rotating-wave approxima¬ 
tion to Eq. (A3). The Markov approximation is based on 
the assumption that within the optical memory time Topt 
all the mechanical dynamics can be neglected except for 
the evolution provided by the free Hamiltonian 
As a result we can write e^^^[xmPm{t—^)] ~ Xm{r)pm{t) 


with Xmir) = e 




^ g other hand. 


the rotating-wave approximation consists in neglecting 
all the terms proportional to and aj^ in the effec¬ 
tive mechanical master equation, under the assumption 
that their rotation at frequency 2Qrn is much larger than 
the rates they are wei ghte d by. After applying these ap¬ 
proximations in Eq. (A3) we are left with an effective 
mechanical master equation given by 


where we have defined the mechanical position quadra¬ 
ture x^ = Om + the photon number in the signal 
mode Ng = trjalasPopt}, and the optical correlation 
function 


s ( t ) = tr{aJase^°p‘'^[a|as/)opt]} - (^-4) 


P^optT 


It is well known that the steady state popt of the DOPO 
is unique (which intuitively comes from the fact that both 
the pump and signal modes have local dynamics leading 
to unique steady states, and the parametric interaction 
preserves that uniqueness), and hence is a relaxing 

map Bam, mapping all optical operators O into the 
steady state, that is, lim^-^oo = tropt{0}popt- 

Thus, the optical correlation function s{r) will always 
decay to zero within some finite memory time which we 
denote by Topt- Hence, in the asymptotic limit we can 
write \imt^ooPm{t - r) = pm{t) = Pm in the 

integral Kernel of Eq. (A3), obtaining an equation for 


pm which is quadratic in the operators Om and therefore 
allows for a Gaussian-state solution [45l |46]. In other 
words, the equations for the first and second steady-state 


Pm — klmPm ’^^mVOMBtg[Xm^ Pm] 

-E Jm^-'^amlPm] + [Pm] 5 

where the heating and cooling rates 7mr± = 
^mVoM^^{^T\'ym=o} coincide precisely with those de¬ 
fined in the main text. This master equation has a very 
simple Gaussian steady state pm corresponding to a dis¬ 
placed thermal state with mean 


{Om) 


iklmVouNg 

iklm ~t' 7m 


VouNg, 


(AS) 


phonon number {6al^6am) = ffm, where is given by 
Eq. (4) in the main text, and {Sa‘^) = 0. We note that, 
starting from a thermal state, the mechanical mode re¬ 
laxes to this steady state with a rate 7eff = 7m (1 + r)^ 
where r = r_— T+is what we have called the cooling 
efficiency in the main text, since the equations of motion 
for the phonon number ffuctuations and the mechanical 
field amplitude are given by 

di{5oJmdamJ) 27m(f ff r) A 27m(^th T—)? 

i^m) ~ [ ^^m 7m (1 ffT)] (^m) H“ ^VOM^m • (A9) 










We have checked that this rate 7 eff is smaller than the de¬ 
cay rate of the optical correlator s{r) for the parameters 
of interest, hence making Markov a valid approximation. 

Let us remar k tha t through out t he work we have been 
using both Eq. (A3) and Eq. (A7) to obtain the steady- 
state moments of the mechanical oscillator. We have 
never observed any notable differences between them, ex¬ 
cept when working extremely close to threshold within 
the semi-classical approach, see the inse t of Eig. 3. In 
these cases, however, the failure of Eq. (A7) can be di¬ 
rectly attributed to the failure of the semi-classical ap¬ 
proach, and not to a failure of the rotating-wave approxi¬ 
mation itself, which indeed is very well satisfied as shown 
by the c-MoP approach. Thus, we conclude that for 
the parameter regime studied in this work (compatible 
with current cWGM resonators) the state of the mechan¬ 
ical oscillator is indeed a displaced thermal state, with a 
phonon number that can only be evaluated once the op¬ 
tical photon number Ng and correlation function s{r) are 
known. 


Appendix B: Semi-classical approach 


column vector 

v { t ) = col ^(aks), (a2), (ak)^ , (B2) 

where the expectation value of an operator A with the 
tilde is defined as {A) = tr{Ae^^^/is}. Taking the deriva¬ 
tive of this vector with respect to r, we find the linear 
system drv{T) = Lv{t), where the matrix M reads 

/ -2 cr ^ \ 

L = 7 o 2a -2(1-fiA) 0 . (B3) 

\ 2cr 0 -2(1 -iA)/ 

It is straightforward to solve this linear system, for ex¬ 
ample by diagonalizing L. We write L = UAU~'^, with a 
similarity matrix U that can be found analytically (but 
its expression is too lengthy to be reported here), and 
a diagonal matrix A containing the eigenvalues of L, 
Ai = — 270 , and A 2,3 = —2^q{1±.\\//S? — cr^). Notice that 
for cr > A the square root becomes imaginary, making 
^2 < 705 and in fact A 2 = 0 at threshold, a = \/l + A^. 
Consequently, we call the region with a > A the critical 
slowing down regime. The solution of the linear system 
is then found as 


The simplest way of obtaining the optical correlator 
s{t) is by using standard linearization on >Copt- In this 
approach, we move to a displaced picture in which the 
large coherent background of the pump mode is removed, 
and then keep terms of the transformed optical Liouvil- 
lian only up to second order in the the bosonic operators. 
The displacement operator D = exp[ep(ap — nJ,)/ 7 o] al¬ 
lows us to move to the new picture, in which the trans¬ 
formed optical state Popt = PoptD evolves according 
to a transformed Liouvillian >Copt = Removing 

terms beyond quadratic order, this transformed Liouvil¬ 
lian can be written as a sum of independent Liouvillians 
for the pump and signal modes, >Copt = Cp Cs^ with 
= 7 o^ap and 


7o-^>C,{.)= [-iAaK+ 2 (Bl) 


with the injection parameter a = Cpxho normalized 
detuning A = As/ 70 . Consequently, the optical steady 
state in the original picture becomes the separable state 
Popt = kp/ 7 o)(ep/ 7 o| c Ps where lep/ 70 ) is a coherent 
state of amplitude 6^/70 and ps is the Caussian state 
satisfying Cs[ps] =0. The latter is completely character¬ 
ized by its first and second moments, which are trivially 
found to be (a^) = 0 , {a\as) = cr^/ 2 (l + A^ — cr^) = TV^, 
and (a^) = cr(l — iA)/ 2 (l + A^ — cr^), where we use the 
usual notation {A) = tr^jAps} for any operator A acting 
on the signal subspace. 

The optical correlation function simplifies to s(t) = 
tr{a|ase^^where we have defined a traceless oper¬ 
ator Ps = (a|as — Ns)ps- Using again the fact that the 
Liouvillian jCg is Caussian, it is simple to evaluate the 
correlation function s{r). To this aim, let us define the 


3 

v{t) = Ue^^U-^v{0) = Y, u, (B4) 

n=l 


where we have defined the initial condition vector 


u= v{0) = col ((a^asa^tts) - (B5) 

(«kks) - {afa,} - (a^yNe), 

and the matrices = 1 / 11 ^ 1 /“^, where ( 11 ^)^/ = dj^din- 
Note that the vector u is formed by fourth order mo¬ 
ments. In order to find them, we simply exploit the 
Caussian structure of £ 5 , which allows us to express mo¬ 
ments of any order as products of moments of first and 
second order. Specifically, concerning third and fourth 
order moments we simply use 

{5a\5a^Q) = = 0 (B 6 a) 

= (5^2) (5a2) + 2(5aka«)^ (B 6 b) 

{5a\5ay) = 2){5a\5as){5ay), (B 6 c) 


where 5as = — (ag). Note that the optical correlation 

function we are looking for is given by the first component 
of the vector, s(r) = [T(r )]i, and the integrals appearing 
in do and d± in Eq. ( |A 6 [ ) can be easily evalu ated due to 
the exponential r-dependence of v{t) in Eq. (B4). 


Appendix C: Cooling by heating via a two-photon 
process 

In the main text we have shown that significant cooling 
can be obtained when working in the resolved sideband 
regime ft = U ^/70 ^ 1 and close to the resonance con¬ 
dition deff = \J AA — a^ = ft 12. We have interpreted such 
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a phenomenon as a “cooling by heating” effect, for which 
we have moved to a new picture defined by the squeezing 
operator S{r) = exp[—ir(a|^ + within the semi- 

classical approach explained above. Here we want to 
explicitly perform an adiabatic elimination of the opti¬ 
cal mode in this picture, what will allow us to get more 
insight into the cooling mechanism. 

Let us first write the master equation in this “squeezed 
picture”. The transformation diagonalizes the Hamilto¬ 
nian in the optical Liouvillian (Bl), turning it into 

7o"'5t(r)£opt5(r)[-] = i[(5effata„ •] (Cl) 

+ (1 + [•] + [•] + iM/Ca, [•] - iM/C^t [•], 


where we have defined the superoperator /Cj['] = 

2J(') J — J^(') — as well as the parameters TVeff = 

(A/^eff — l)/2 and M = cr/2(5eff. Note that the JC 
terms rotate at frequency 2Jeff 5 and hence are highly 
suppressed when we work within the cooling condition 
Jeff ~ ^/2 and with M/2 Jeff = cr/dJ^j^ <C 1 (rotating- 
wave approximation). Therefore, as mentioned in the 
main text, in this picture the signal field is turned into 
a bosonic mode with oscillation frequency Jeff, and at 
thermal equilibrium with occupation TVeff- On the other 
hand, the photon-number operator is transformed into 
5t(r)ata,5(r) = N,s + (2iVeff + l)ata, + iM{al - aff), 
and hence the optomechanical interaction can be approx¬ 
imated by 


S\r)alasS{r){am + - 4^am)- (C2) 


within the rotating-wave approximation as long as 
7/OM(2A^eff + 1) = VOM^/^eff 1- 

Hence, within these conditions, the transformed state 
p = S'^{r)pS{r) evolves according to a master equation 
that we write as 


dtp — Cs[p\ + ^m[p\ + ^om[p\ (03) 

with 

^AP] = [-i 7 o^^effala 5 ,p] (C4a) 

+7o(Aeff + l)^as [p] + 70 ^eff^Qt [p], 

~ (C4b) 

H“7m('^th H“ f)^arn,[p] T 7m^thL^Ql^ [p]: 
ComIp] = [VLmPoMM{ala\^ - a^^a^),p], (C4c) 

where TVeff and M are defined in the main text. The 


structure of thi s ma ster equation is similar to the orig¬ 
inal one, Eq. with the only difference that the 

optical Liouvillian is replaced by £ 5 , corresponding to a 
single-mode at finite temperature, and the optomechani¬ 
cal interaction a|as(a^ + aj^) by iM(a^aJ^ — a|^a^). The 
adiabatic elimination of the optical mode can be carried 
out exactly in the same way as we did in Appendix [A| 
and under the cooling condition Jeff = f 2/2 it would lead 
to the heating and cooling rates 


r_ ^CMhiialal'^ps}, 

(C5a) 

r+ ^CMhv{al^alps}, 

(C5b) 


where C = ^^PoM/ 7 rn 7 o is the bare optomechanical co- 
operativity, and ps is a thermal state with mean photon 
number TVeff- The cooling efficiency is then given by 

r = r_-r+ = ^CMhT{[alal^]p,} = 2 CM^(N,f,+ ^. 

(C6) 

The “cooling by heating” effect is clearly seen because 
the cooling efficiency increases with the effective ther¬ 
mal photon number TVeff- But it is important to note 
that this enhancement of the cooling efficiency is a di¬ 
rect consequence of the commutator appearing in the 
trace, contributing as [a^, = 4a|as + 2 , which in turn 

comes from the fact that the effective optomechanical in¬ 
teraction i(a^aj^ — a|^a^) corresponds to the exchange 
of phonons with pairs of photons. In the usual sideband 
laser cooling scenario, the effective optomechanical in¬ 
teraction is bilinear, e.g. i(asaj^ — a|a^), meaning that 
the commutator in the expression above is replaced by 
[ug, a|] = 1 , and hence the thermal photonic background 
does not enter the cooling efficiency. 

Let us finally note that the fundamental limit can be 
written as 

- _ T^ _ tr{a| UgPg} _ -^eff Aeff>i -Neff 

T tr{[a 2 ,a|^]ps} 27Veff +1 2 

(C7) 

which increases linearly with the effective thermal pho¬ 
ton number. Hence, as explained in the main text, the 
“cooling by heating” mechanism is optimized by finding 
a proper trade off between the increase in the cooling 
efficiency (“good noise”) and the increase in the fun¬ 
damental limit (“bad noise”). It is to be noted that 
within the usual sideband laser cooling, any thermal 
background will still contribute to this fundamental limit, 

ffpL = tT{alasPs}/tr{[as, al]ps} Neff, but, as ex¬ 

plained above, it provides no enhancement of the cooling 
efficiency T. In other words, in standard sideband cooling 
the thermal background acts only as “bad” noise. 

We emphasize that these expressions for T and npL 
agree with the ones provided in the main text, which have 
been first calculated exactly within the semi-classical ap¬ 
proach, and then approximated to leading order in 1/Q‘^ 
for Jeff = n/2. 


Appendix D: c-MoP approach to the optical problem 
and scalings at the critical point 

Despite the analytical and physical insight that it pro¬ 
vides, the semi-classical approach suffers from several is¬ 
sues, most importantly its divergent character at thresh¬ 
old, which shows that it cannot be trusted when working 
close to such point. Unfortunately, there is no systematic 
way of checking within the formalism itself where exactly 
it fails. This question can only be answered by compar¬ 
ing it to a more accurate approach. To this aim, we 
have applied our recently developed self-consistent Mori 
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FIG. 4. (Color online) (a) Steady state photon number Ns = tr{aJasPopt} at the critical point as a function of the normalized 
detuning A. (b) Absolute value of the normalized optical correlation function s(t) at the critical point for different values of 
detuning as of function of the normalized time 7optT, with 7opt = 7 o^dc (1 + A). In both figures the parameters are 70 = 1 and 
?7 dc = O.Ol, and we have obtained them by applying Gaussian c-MoP theory to the optical problem. 


projector (c-MoP) theory [24l|25]. In particular this ap¬ 
proach has allowed us to characterize optical steady-state 
observables such as the photon number Ng or the corre¬ 
lation function s{r) in all relevant parameter space, in¬ 
cluding threshold. 

We detailed the application of c-MoP theory [24] to 
the DOPO problem in m, including its combination 
with the Gaussian-state approximation that we use in 
this work, which was shown to be quite accurate both for 
steady-state quantities and dynamics. Let us now briefly 
introduce such an approach here for completeness, keep¬ 
ing in mind that details can be looked up in [25]. Our 
main goal consists in finding the optical correlation func¬ 
tion (§, for which we need to solve the dynamics gener¬ 
ated by the optical Liouvillian >Copt, that is, the DOPO 
dynamics 


Popt — ^^optPopt- 


(Dl) 


This is a nonlinear two-mode problem with no analytic 
solution and whose direct numerical simulation becomes 
unfeasible for moderate photon numbers already. In 
contrast, c-MoP theory works with the reduced states 
of the pump and signal modes, pp{t) = tr^jpopt} and 
ps{t) = trpjpopt}? respectively, whose coupled dynam¬ 
ics are approximated by the set of nonlinearly coupled 
equations 


m + I [at2(a,)(t) - a2(at)(t),p,(t)] 

+ (|)'{[a^/i,(i)]+H.c.}, (D2a) 

Pp{t) = CpPp{t) + I [ap{al^){t) - al{al){t),pp{t)] 








n=l 


+ H.c. } , (D2b) 


hs{t) = {--fp N Cs)hs{t) N JCs{t,t)ps{t), (D2c) 

hpA'^) = + Cp)hp^n{t) + t)pp{t), (D2d) 


where hs{t) and {/?'p,n(t)}n=i, 2,3 are auxilliary operators 
acting on the signal and pump subspaces (introduced to 
turn the c-MoP equations into ordinary differential equa¬ 
tions, since originally they have an integro-differential 
structure [25]), and we refer to [25] for the definitions of 
the superoperators £ 5 , Cp^ /Cs(t,t'), and /Cp,n(L^ 0 - 
noting by Dg and Dp the dimensions of the truncated 
Hilbert spaces of signal and pump in a num erical simula¬ 
tion, we see that the original problem ( jPlj ) requires solv¬ 
ing 8i Dg X Dp system, whereas there are only 2Dg + ADp 
c-MoP equations (or even less by exploiting the quadratic 
or Gaussian structure of the pump equations). 


Nevertheless, even though c-MoP allows to get numer¬ 
ical insight in a larger region of parameter space, the 
simulation of problems with very large photon numbers 
(such as the ones we work with close to threshold) is still 
challenging. It is in those regions where a Gaussian-state 
approximation becomes extremely useful. As the name 
suggests, such approximation consists in assuming that 
the reduced signal and pump states are Gaussian, mean¬ 
ing that they are completely characterized by first and 
second order moments. Under such circumstances, we 
can a pprox imate third and fourth order moments of pg 
as in (B 6 ), and the c-MoP equations provide a closed 
set of nonlinear equations for the first and second or¬ 
der moments of the operators Ps(t), Pp{t), hg{t) and 
{^p,n(^)}n=i, 2 , 3 - The stcady-statc moments can then be 
efficiently found simply by finding the stationary solu¬ 
tions of these equations. 


As an example, in Fig. if we show the steady-state 
photon number Ng at the critical point {x = 1 ) as a 
function of the normalized detuning A. It shows a clear 
linear dependence for A > 1 which, together with the 
well-known scaling with the down-conversion cou¬ 
pling usumli], provides an overall Ng oc (1 + A)/7^dc 
scaling of the signal photon number at threshold. The 
knowledge of this scaling plays an important role for the 
determination of the conditions under which mechanical 
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backaction on the optics can be neglected, as we have 
explained in Sec. lYg 

Let us now explain how the optical correlation function 
can be evaluated within this framework. First, note that 
we can rewrite it as 


s(t) = Ns (tr{a|a5Z/(T)} - Ns 


(D3) 


where u{t) = can be interpreted as an oper¬ 

ator with evolution equation z> = £opt^ and initial con¬ 
dition z^(0) = alaspopt/Ns. Since this evolution equa¬ 
tion is fo rmally equivalent to the optical master equa¬ 
tion ( |D1| ), we can apply c-MoP theory directly on z^(t), 
approximating it by a separable operator Usit) 0 
with = tTs{iy{t)} and z/s(t) = tTp{u{t)} evolving 

according to equations ( |D2| ) with pj replaced by Uj . Un¬ 
der a Gaussian approximation for Usit) similar to (B6) 
but with expectation values defined with respect to Ug (t), 


the evolution equations for the first and second moments 
of Up{t), hs{t), {hp,n(^)}n=i,2,3, and their Hermi- 

tian conjugates (note that v is not Hermitian) form a 
closed nonlinear system which we can solve again effi¬ 
ciently. Note that the initial conditions for these mo¬ 
ments, e.g., tr{a|asZ/(0)} = tr{a|asa|asPopt}/^s, are 
found from the Gaussian c-MoP steady-state solutions 
as explained above. 

In Fig. i) we show the evolution of the absolute value 
of the correlation function s{t) at the critical point and 
for different values of the normalized detuning A. Time is 
normalized to [7 o?7dc(1 + and hence the fact that 

all the curves decay on the same time scale proves that 
the optical relaxation time scales as 7opt = 7 o7dc(1 + 
at threshold. This again plays a fundamental role when 
proving that mechanical backaction is negligible, as we 
discussed in Sec. I VII 
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